Continuum QRPA in the coordinate space representation 



K. Hagino 1 and H. Sagawa 2 
Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 
2 Center for Mathematical Sciences, University of Aizu, Ikki-machi, 
Aizu-Wakamatsu, Fukushima 965-8580, Japan 



Abstract 
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F ^ \ We formulate a quasi-particle random phase approximation (QRPA) in 

P_H \ the coordinate space representation. This model is a natural extension of the 

t*"*- ' RPA model of Shlomo and Bertsch to open-shell nuclei in order to take into 

account pairing correlations together with the coupling to the continuum. We 
apply it to the 120 Sn nucleus and show that low-lying excitation modes are 
significantly influenced by the pairing effects although the effects are marginal 
■^-j- \ in the giant resonance region. The dependence of the pairing effect on the 

parity of low-lying collective mode is also discussed. 
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I. INTRODUCTION 



The random phase approximation (RPA) has provided a convenient and useful method 
to describe excited states of many-fermion systems. There are a number of ways to formu- 
late the RPA |l]-|9| . In practical point of view, we particularly mention here a configuration 
space formalism and a response function formalism. The configuration space formalism di- 
agonalizes a non-hermitian matrix denoted often by A and B matrices which are constructed 
in the model space of 1 particle- 1 hole (lp-lh) states. In contrast, the response function 
formalism is based on the linear response theory and solves a Bethe-Salpeter equation for the 
response function, often in the coordinate space The response function formalism 

can be formulated in various representations of basis states, and the connection between 
these two methods can be made by expressing response functions in the configuration space 

i- 

The response function formalism of RPA becomes particularly simple when the interac- 
tion is a zero-range contact force. In that case, the excitations to particle continuum states 
can be treated exactly by solving the single-particle Green function in the coordinate space. 
This method was first developed by Shlomo and Bertsch [|J and subsequently applied to 
self-consistent calculations of nuclear giant resonances with Skyrme interaction by Liu and 
Van Giai [011. Recently it has been extensively used by Hamamoto, Sagawa, and Zhang to 



discuss giant resonances of neutron-rich nuclei, where the continuum effects play an essen- 
tial role due to a much lower threshold energy compared with /5-stable nuclei ]TT],[n^]. An 
extension to the three dimensional space has also been carried out recently by Nakatsukasa 



and Yabana for the study of atomic clusters fl3fl . The coupling between particle-particle 
continuum was also studied in 11 Li by Bertsch and Esbensen [H,|T5||. 

Although it has been well known that the pairing correlation plays an important role in 
the ground state of open-shell nuclei, it has been neglected in applying the response function 
formalism to describe excited states of atomic nuclei until very recently ||16|| . So far a simple 
filling approximation has been employed in open-shell nuclei in order to simulate the pairing 
correlations. However, it is important to take into account the pairing effects consistently 
for the study of excited states of open-shell nuclei, and thus the quasi-particle RPA (QRPA) 



should be used instead of the RPA fL7 |. 



The aim of this paper is to generalize the formalism of Shlomo and Bertsch to the QRPA 
and discuss the effects of pairing correlations on excited states of open-shell nuclei. In this 
paper, we study only cases where the BCS approximation works well and thus all states in 
the pairing active space are bound, leaving out as a future work a self-consistent treatment 
of the pairing effects in the continuum using the Hartree-Fock-Bogoliubov + RPA theory 
P~8|-PD|| . Our method is closely related to that in Ref. |TBL although details are somewhat 
different. The paper is organized as follows. In Sec. II, we briefly review the formalism of 
Shlomo and Bertsch and extend it to the QRPA. In Sec. Ill, we apply the formalism to 
the isoscalar (IS) monopole, quadrupole, and octupole modes as well as the isovector (IV) 
dipole mode of excitation of the 120 Sn nucleus and discuss the effects of pairing correlations 
on the excited states. A summary is given in Sec. IV. 
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II. LINEAR RESPONSE THEORY FOR OPEN-SHELL NUCLEI 



We begin with the configuration space formalism of the RPA theory and then make a 
connection to the response function method. The RPA equation can be given in a compact 
form 



A B \ fX\ (X 

UJ 



B* -A* \Y \Y 



(1) 



where X and Y are the forward and the backward amplitudes, respectively. For a residual 
interaction of delta-type contact force, v res (r 1 ,r 2 ) = v((ri + r 2 )/2) 5(r*i — r 2 ), the matrix 
elements of A and B for the L-multipole mode read 

A ph ,p'h' - (e p - eh)$ P h,p'h> = B phjp/h/ = I(php'ti) (j p l p \\Y L \ \j h l h ) (j p il p , \\Y L \\j h il h i) ^ - , (2) 

where p(h) denotes a particle (hole) state and \Y L \ \j'V) is the reduced matrix element. 
The radial integral I is given by 

/■OO rhr 

I(php'h') = J o ^v(r)<f> p (r)<t> h (r)Mr)Mr), ( 3 ) 

where <f>(r) is a single-particle radial wave function. We have absorbed the overall factor 
(— 1) L in front of the B matrix by redefining the sign of the backward Y amplitudes. 
The integral (^) can be computed by discretizing it as 

Ar 

I{php'ti) w — v ( r k)4>p(rk)4>h(r k )(f) p >(r k )(t) h/ (r k ), (4) 
k r k 

where Ar is the spacing of the radial coordinate. Since the interaction is then given as a 
sum of separable form, the RPA frequencies uj can be obtained by solving the generalized 
RPA dispersion relation det(l — IIo(a;)x) = 0, where the matrices Hq(u) and x are given 
by 



no(*\j;w) = -J2D ph (i)D ph (j) 

vh \ 



1 1 

jih v . e h — L}-ir) e p - e h + uj - ir] / 



(5) 



Ar 

X(hj) = X(i) kj = — v ( r i)i (6) 

respectively, in the coordinate space representation. Here, r] is a infinitesimal real number, 
and D p h is given by 

D ph {i) = (j) P {ri)4>h{ri){j p l p \\Y L \\j h l^ -^== • (7) 

Note that n (cc;) is nothing but the unperturbed response function in the linear response 
theory. Here we introduce the RPA response function which obeys a Bethe-Salpeter equation 

Rrpa = n + TIoXRrpa- (8) 
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This equation can be solved in the coordinate space by matrix inversion as 

rw(z,,7» = £(i - u (u)x)Tp (Kj;^). (9) 

k 

The response of the system to an external field V ext (r) = V ex t{r)Yi,M{r) is then given by |7| 
S(u) = J2\(f\V ext \0)\ 2 5(E f -Eq-u) = hm J dr t J drjV^r^URPAihj^V^Tj). 

(10) 

The exact treatment of the continuum effects can be achieved by eliminating the sum 
of the particle states in the free response function (|J) using the complete set of the wave 
function HI. This leads to 



h 



i 



jplp 



2L + 1 



x n 



h-e h 



+ 



1 



uj 



17] h — €h + to — if] 



(11) 



where h is the single-particle Hamiltonian, and the single-particle Green function is given 
by 



1 



h — eh ± uj — it] 



2m u(r<)w(r>) 



h 1 



W 



(12) 



Here, u and w are the regular and irregular solutions of the Hamiltonian h at energy ^fuj, 
and W is the Wronskian given by W = uw' — wu' . 

In the application of the formalism to nuclear systems, we use the proton-neutron for- 
malism which can properly take into account the couplings between the IS and IV modes of 
excitation j!!]. The Bethe-Salpeter equation given by Eq. (P) is then generalized to be 
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(13) 



For simplicity of the notation, we do not present below the formalism in the two-component 
form, although we indeed use it in the actual calculations given in the next section. 

In the presence of paring correlations, the elementary excitation is a two quasi-particle 
excitation, rather than a particle-hole excitation. A generalization of the RPA to the QRPA 
is given in Ref. together with relevant A and B matrices in the QRPA equation. A free 
response function in the QRPA in the configuration space representation can be constructed 
in a similar way to the RPA, 



a</3 



1 



+ 



E a + Ep — uj — if] E a + Ep + uj — if] / 



(14) 



with 
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D a p(i) = <j>a{ri)4>p{ri)(jJ a \\Y L \\jplp) 



UgVp + (-; rv a up 
V2L + 1 ' 



1 + 5, 



a,0) 



-1/2 



(15) 



where u 2 is the BCS occupation probability and u 2 a 
energy given by E a = 



1 — v\. E a is the quasi-particle 

(e a — A) 2 + A 2 , where A and A Q are the chemical potential and the 
pairing gap, respectively. In the BCS approximation, <p a is an eigenfunction of the single- 
particle Hamiltonian h with an eigen-energy e a . Since the quasi-particle energy E a is not 
an eigenvalue of h in general, it is not straightforward to introduce the single-particle Green 
function ( |T2D in the QRPA free response function. However, when the pairing gap is zero 
for states k outside the pairing active space, the quasiparticle energy becomes an eigenvalue 
of the single-particle Hamiltonian, i.e. E^ = e k — A, i>& = 0, and Uk=l- We therefore consider 
separately excitations among states within the pairing active space and those from the inside 
to the outside of the active space. To the latter model space, we apply the same procedure 
as the RPA response function. The free response function in the BCS approximation (14) 
thus becomes 



a<3 V 



J2(jJa\\Y L \\j k lk) 



E a + E/3 — uj 
1 



+ 
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it] E a + Ep + uj - irj 
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+ 



2L + 1 
1 



E a + h — A — uj — ir\ E a 



h — A + uj — it] 
1 



+ 



E a + ep — A — uj — ir\ E a + ep — X + uj — irj / 



(16) 



where the summations of a and (3 are restricted to the states within the pairing active 
space. The last term in Eq. ([IB] ) is a correction for a double- counting of excitations within 
the pairing active space, which stems from the substitution of the completeness relation 
Ylk\4>k) i&kl — 1 — Hb\4 > p){4 > p\ m Eq. QI4]). Note that, without the pairing correlations, 
E a = X — e a and v a = 1 for e a < A, and the BCS free response function fllB]) is identical 
to the corresponding response function ( ]TT| ) in the RPA. With this unperturbed response 
function (|T6|), the QRPA response function is obtained again by solving the Bethe-Salpeter 
equation (B) as in the RPA theory. 



III. CONTINUUM QRPA EXCITATIONS IN 120 SN 



We now apply the continuum QRPA formalism to 120 Sn and make a comparison between 
the QRPA and the RPA. We choose this system because it is a typical open-shell nucleus and 
also because it is a sub-shell closure nucleus where the RPA can be applied unambiguously 
without using the filling approximation for valence nucleons. The single-particle wave func- 
tions <p and the single-particle energies e are obtained by solving the Shrodinger equation 
with a Woods-Saxon mean-field potential. As a residual interaction v res , we use the to and 
£3 parts of the Skyrme residual interaction, which is obtained from the second derivative of 
the Skyrme energy functional with respect to the proton and the neutron densities. The 
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ground state density to be used in the density-dependent t 3 part of the residual interaction 
is generated from the single-particle wave functions <j). We use the same parameters as those 
in Ref. p| for the mean-field potentials and the residual interaction, i.e., to = —1100 MeV 
fm 3 , t 3 =16000 MeV fm 6 , xo=0.5, £3=1, and 7=1 in the standard notation of the Skyrme 
functional pi| . Since this model is not self-consistent, we renormalize v res for the QRPA 
and the RPA calculations, respectively, so that the spurious IS dipole mode appears at zero 
excitation energy. The resultant renormalization factors are k = 0.638 and 0.660 for RPA 
and QRPA, respectively. For the pairing of neutrons, we use a schematic state-independent 
constant pairing gap A Q = A = 1.392 MeV, which is estimated from the experimental 
binding energies of neighboring nuclei. The pairing active space which we adopt includes 
the levels up to the A = 50 major shell as well as the lg7/2, 2ds/2, 2d3/2, 3si/2, and lhn/ 2 
levels. The single-particle energy for the valence levels are shown in Table 1, together with 
the BCS occupation probabilities. The proton pairing gap is set to be zero due to the magic 
number Z = 50. 

Figs. 1-4 show the strength function S(u) for the isoscalar (IS) monopole, IS quadrupole, 
IS octupole, and the isovector (IV) dipole modes of excitations, respectively. The external 
field V ex t(r) for each modes is given by r 2 , r 2 , r 3 , and rr 2 , respectively. The lower panels 
show results of the QRPA calculations, while the upper panels show results of the RPA 
calculations without the pairing correlations as a comparison. For the latter, the levels up 
to 3si/ 2 state are occupied with the occupation probability v 2 of 1, while those above lhn/2 
state are unoccupied with v 2 = (see Table 1). Arrows indicate the positions of (Q)RPA 
states below the threshold, where the strength function has no width. 

For the monopole excitation shown in Fig.l, the RPA does not show any low-lying mode, 
although the experimental second lowest + state in the 120 Sn nucleus is observed at the 
excitation energy 1.874 MeV, which is attributed to a paring vibration mode. The second 
lowest + state in the QRPA is found at rather low energy 2.9 MeV, and may couple to 
the pairing vibrational mode. The IS giant monopole states are seen in the energy region 
17—28 MeV in both the QRPA and RPA calculations. The transition densities for the IS 
monopole mode at two different energies, i.e., u = 2.9 MeV and 21.6 MeV, are shown in 
Fig. 5. We find that the former shows a characteristic behaviour of the pairing vibration 
mode while the latter shows the compressional character. For the quadrupole mode shown 
in Fig. 2, the lowest RPA state is at 5.2 MeV, in comparison with the experimental value of 
1.17 MeV. The transition strength is B(E2:0 + -> 2+)= 0.031 e 2 6 2 in the RPA calculation, 
while the experimental value is 0.2 e 2 6 2 . If one takes the pairing correlations into account, 
the low-lying 2 + RPA states goes substantially down and the lowest QRPA state appears 
at 2.3 MeV having the transition strength B(E2:0 + — > 2 + )= 0.107 e 2 b 2 , which is much 
closer to the experimental value. In both cases, the main peak of the IS giant quadrupole 
resonance (GQR) appears at around 13 MeV exhausting most of the sum rule strength. 



The experimental GQR is also observed at the excitation energy around 13MeV [22 . As for 



the octupole mode of excitations in Fig. 3, the experimental value of the lowest 3 state 
is at 2.4 MeV, while the RPA and the QRPA lead to 1.4 MeV and 3.0 MeV, respectively. 
The experimental value for the corresponding transition rate is B(E3:0 + — > 3~) = 0.09 e 2 b 3 , 
while it is 0.159 e 2 6 3 and 0.0771 e 2 b 3 in the RPA and the QRPA, respectively. The RPA 
without pairing correlations underestimates the lowest 3~ state energy and the QRPA is 
more satisfactory compared with the experimental value. We summarize in Table 2 the 
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results of the RPA and the QRPA calculations together with the experimental data for the 
lowest 2 + and 3~ states. Giant octupole resonances (GOR) are observed in Fig. 3 at the 
excitation energy 23 MeV in both cases. Results for the IV giant dipole resonance (GDR) 
are shown for RPA and QRPA in Fig. 4. As is the same as in Figs. 1-3, the structure of 
GDR is not much disturbed by the effect of the pairing. 

Notice that the effect of the pairing on the low-lying quadrupole mode is opposite to that 
on the octupole mode. For the former, the energy of the lowest 2 + state becomes smaller 
due to the pairing, with enhancement of the B(E2) value. On the other hand, the QRPA 
raises the energy of the lowest 3~ state and the B(E3) value is hindered. In both cases, 
the QRPA gives better agreement with the experimental values compared with the RPA. 
These features can be understood as follows. For excitations with even parity, such as the 
quadrupole excitations, transitions between the same levels are allowed in the presence of the 
pairing correlations. This enhances the collectivity of low-lying states, resulting in a smaller 
excitation energy and a larger transition strength. On the contrary, such excitations are 
not allowed for odd parity excitations, e.g., the octupole mode. In that case, the dominant 
excitation is from one major shell to another regardless of the pairing correlation. Since the 
particle-hole excitations are weakened by the factor v 2 in the presence of the pairing, the 
QRPA lowers the collectivity of low-lying collective excitations with odd parity. 

In contrast to the low-lying modes of excitations, in general, high-lying modes are much 
less sensitive to the paring correlations. This is the case for all the modes of excitations 
shown in Figs. 1-4. This fact can be seen more transparently if one smears the strength 
functions with a finite width. We show in Figs. 6-9 the RPA and QRPA results with a finite 
value of 7] =0.5MeV in Eq. (|16|). As one clearly sees, the two results resemble each other in 
the giant resonance region at energies above the threshold near 10 MeV. These results are 
natural consequences of the pairing correlations since the configurations outside the pairing 
active space are the main p-h configurations for the giant resonances and thus the pairing 
effects should play a minor role. 

A conventional approximation to treat the continuum effect is to put a nucleus in a 
box and discretize the continuum states by imposing a boundary condition at the edge 
of the box. A model space of these states is usually truncated at the maximum single- 
particle energy e max . Since our continuum QRPA method can treat the coupling to the 
continuum exactly, it is interesting to compare our method with the box discretization 
approximation. Fig. 10 shows a convergence property of the box discretization method for 
the IS quadrupole response. We take the box size of R max =10 fm, and smear the strength 
function with 77=0. 5 MeV for the representation purpose. We have checked that the results 
do not change significantly even if we use a larger value of R ma x- The solid line is the result 
of the continuum QRPA method with the exact treatment of the continuum effect, which is 
the same as in Fig. 6. The dot-dashed, dashed, and thin solid lines are results of the box 
discretization method with truncation energy at e max = 10, 50, and 100 MeV, respectively. 
In the upper panel, we use the same residual interaction for all the calculations. As one can 
see, the convergence with respect to the maximum energy of continuum states is extremely 
slow: even with e max = 100 MeV the peak positions are not reproduced. The peak height 
for the lowest energy peak is not reproduced, either. We note that the truncation of the 
model space shfts the position of the spurious IS dipole state. It appears at 6.0, 4.4, and 
3.7 MeV for e max = 10, 50, and 100 MeV, respectively. In the lower panel, the renormalized 
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residual interaction is used for each value of e max so that the IS dipole mode appears at zero 
energy. The renormalization factors are k= 0.968, 0.776, and 0.732 for e max = 10, 50, and 
100 MeV, respectively. This procedure drastically improves the convergence. In the case of 
the truncation energy e max = 50 MeV (dashed line), the exact solution is almost reproduced. 
The result converges at this energy and the shape of the strength function is not altered 
even when the states up to e max = 100 MeV are included. This study clearly indicates the 
importance of the self-consistency between the space truncation and the renormalization of 
the residual interaction. If this self-consistency is imposed, the box discretization method 
works well. It should be emphasized, however, that the escape width of resonances can be 
obtained only when the single-particle continuum states are treated exactly as we do in this 
article. 



IV. SUMMARY AND DISCUSSIONS 

We proposed a QRPA model in the coordinate space to take into account the continuum 
effects. It is a generalization of the formalism of Shlomo and Bertsch for the continuum RPA 
to the QRPA. We treated separately the p-h excitations within the pairing active space and 
those between the active space and the non-active space. For the former, we explicitly used 
the two quasiparticle configurations in the response function in the coordinate space, while 
for the latter we use the single-particle Green function in the coordinate space representation 
taking into account the coupling to the particle continuum properly. We applied the formal- 
ism to 120 Sn and showed that the pairing correlations enhance the collectivity of the positive 
parity low-lying states, while that of the negative parity low-lying states are hindered. We 
found that the QRPA is more satisfactory in reproducing the experimental data of the en- 
ergy of low-lying modes of excitation, while the giant resonances are not much affected by 
the pairing correlations. We would like to point out that these results however do not jus- 
tify the RPA model without the pairing correlations for the study of giant resonances in 
open-shell nuclei, especially near the drip-lines. Since the single-particle energies and wave 
functions could be different in principle due to the density dependence of the mean-field, 
the response properties even in the giant resonance region might be affected by the pairing 
correlations. We therefore advocate to use the QRPA all through the excitation region in 
open-shell nuclei. 

There are many possible applications of the continuum QRPA method presented in 
this paper. One of them is a neutrino-nucleus scattering, e.g. 12 C(z/, u') l2 C* |2B| and 
12 C(z/ e , e~) 12 N 1 24}j , where the RPA and the QRPA have been one of the standard meth- 
ods for theoretical investigations. Another application will be the excitations of drip-line 
nuclei. Because of a low energy threshold, the responses would be very sensitive to the pair- 
ing correlations. In these cases, the pairing active space includes particle continuum states 
and the BCS approximation should be carefully examined. Also, the residual interaction 
in the particle-particle channel, which we neglected in this paper, might have to be taken 
into account. A more sophisticated treatment of the pairing interaction is provided by the 
Hartree-Fock-Bogoliubov (HFB) theory, which consistently describes couplings between the 
particle-hole and the particle-particle channels. It would be an interesting future work to 
develop a continuum QRPA theory based on the HFB approximation. 
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TABLES 

TABLE I. Neutron single-particle levels near the Fermi surface for 120 Sn obtained with a 
Woods-Saxon potential. The occupation probabilities v 2 are calculated in the BCS approxima- 
tion with a constant pairing gap A = 1.392 MeV. The pairing active space includes those levels 
shown in this table as well as the levels up to the N = 50 major shell. The neutron Fermi energy 
A n is -8.149 MeV. 



level energy (MeV) occupation probability 
lg 7/2 -11.184 0.954 
2d 5/2 -11.145 0.953 
2d 3/2 -9.280 0.815 
3s 1/2 -9.048 0.771 
lh 11/2 -6.949 0.173 



TABLE II. Comparison of the RPA and the QRPA calculations with the experimental data 
for the lowest-lying 2 + and 3~ modes of excitation of 120 Sn. 





E 2 + (MeV) 


B(E2)(e 2 b 2 ) 


E 3 _ (MeV) 


B(E3)(e 2 b 3 ) 


RPA 


5.2 


0.031 


1.4 


0.159 


QRPA 


2.3 


0.105 


3.0 


0.0771 


Expt 


1.17 


0.2 


2.4 


0.09 
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FIG. 1. Comparison between the RPA (the upper panel) and the QRPA (the lower panel) for 
the isoscalar monopole mode in the 120 Sn nucleus. Arrows indicate the position of QRPA states 
below the threshold. 
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FIG. 2. Same as fig. 1, but for the isoscalar quadrupole mode. 



30 



13 



1 i 1 r 

IS Octupole 
RPA 



l YY i W 





QRPA 







■Y I Y w vi 



1 



1 








25 



5 10 15 20 

CO (MeV) 

FIG. 3. Same as fig. 1, but for the isoscalar octupole mode. 
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FIG. 4. Same as fig. 1, but for the isovector dipole mode. 
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FIG. 5. Transition densities for the isoscalar monopole mode obtained in the QRPA (the solid 
line). The upper and the lower panels are for uj = 2.9 and 21.6 MeV, respectively. The neutron 
contribution is plotted separately by the dashed line. 
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FIG. 7. Same as fig. 6, but for the isoscalar quadrupole mode. 
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FIG. 8. Same as fig. 6, but for the isoscalar octupole mode. 
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FIG. 10. Convergence property of the box discretization method for the continuum coupling 
for the isoscalar quadrupole mode of 120 Sn. The box size for the discretization is taken to be 10 fm, 
and results are smeard with 77=0.5 MeV. The thick solid line is the result of the continuum QRPA 
calculations which treat the continuum effects exactly. The upper and lower panels are obtained 
without and with imposing the self-consistency condition of the residual interaction, respectively. 
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